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Abstract 

In this paper the scalar-tensor theory of gravity is assumed to describe the 
evolution of the universe and the gravitational scalar (j) is ascribed to play 
the role of inflaton. The theory is characterized by the specified coupling 
function uj{4>) and the cosmological function X((p). The function \(<fi) is 
nearly constant for < cp < 0.1 and A(l) = 0. The functions \(<f>) and u>((/)) 
provide a double- well potential for the motion of <fi{t). Inflation commences 
and ends naturally by the dynamics of the scalar field. The energy density 
of matter increases steadily during inflation. When the constant V in the 
action is determined by the present matter density, the temperature at 
the end of inflation is of the order of 10 1 GeV in no need of reheating. 
Furthermore, the gravitational scalar is just the cold dark matter that men 
seek for. 

PACS: 04.20.-q; 98.80.Cq. 
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I. INTRODUCTION 



The idea of inflation PJ, as well known, gives not only the solution for the horizon 
problem and flatness problem that plagued the standard big bang model, but also the 
explanation of the formation of large-scale structure of the universe from the evolution of 
the primordial density perturbations. According to current approach, during the evolu- 
tion of the very early universe, one usually assumes that the energy density happens to 
be dominated by some vacuum energy and then comoving scales grow exponentially or 
quasi-exponentially, and the vacuum energy driving inflation is generally assumed to be 
associated to some scalar field, the inflaton, which is displaced from the minimum of its 
potential. There results a dilemma. The level of density and temperature fluctuations ob- 
served in the present universe, Sp/p ~ 10~ 5 , require the inflaton potential to be extremely 
flat, and this is in contrast with the requirement that the coupling of the inflaton field 
cannot be too small otherwise the reheating process, which converts the vacuum energy 
into radiation at the end of inflation, takes place too slowly and insufficiently |2| 

The futile effort to identify the inflaton in particle sector and the fact that gravitational 
force is the only long-range force governing the evolution of the universe lead us to consider 
the inflation as a pure gravity effect and the gravitational scalar in the scalar-tensor gravity 
of Bergmann, Nordtvedt and Wagoner || as the inflaton field. 

In section |J we se ^ down the basic formalism. In the action of the scalar-tensor 
gravity theory, there are a coupling function u)(<f>) and a cosmo logical function A(0). The 
function X(4>) is nearly constant for < < 0.1 and A(l) = 0. The functions A(0) and 
uj{4>) provide a double- well potential for the motion of <f>(t). In the Lagrangian there are 
three constants £, (3 and Y in u ((/)), A(0) and the coupling term respectively. The £ and 
(3 are determined by the large-scale structure data, while V is determined by the matter 
content of the universe, i.e., the temperature of CMB and the density of non-relativistic 
matter together. We shall have a model with (£,/?, T) = (7.5, 1.7 x 1(T 16 , 0.06362). The 
basic cosmological equations for the Robert son- Walker metric are given in section |T|. 



In section |TV| we discuss the very early universe. The universe is created from nothing 
but gravity. After a quantum era, we set the initial condition at the Planck time. The 
inflation commences naturally because of the pressure of the scalar field <fi is negative 
and equals in magnitude with the density p^ of (ft. The gravitational scalar field is just 
the inflaton field. 
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In section [V[ we discuss the cosmological solutions after inflation. Section [V| proposes 
a realist model. Section VII discusses the primordial nucleosynthesis. 



II. THE SCALAR-TENSOR GRAVITY THEORY 

We have proposed such a scalar-tensor gravity theory ||. It is characterized by an 
action: 



A = J d 4 x^—g 



\2 
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in which L m is the Lagrangian of matter, • = 4> tlJ (p ,cr , u • 4> = w (U 0' M , m m is the four- 
velocity, ,At = <90/<9x M , and T is a constant of mass dimension 0. The T— term describes 
the coupling of the field with matter which is described as an ideal fluid. The coupling 
function u(4>) and the cosmological function A(0) are given as 

2^(0) + 3 = -^-, (2) 

1 - 

A(0) = 2£9(l-0-0Zn0), (3) 
where £ and (3 are two dimensionless constants. In a model, we can set (3 = 1.7 x 10~ 16 



and £ = 7.5 by the data of COBE DMR observations (we shall see in section |Vl|). This 
value of (3 is extremely small indeed, yet there is no fine-tuning problem, since it is not 
a parameter in gauge theory but a number to characterize the cosmological function A. 
We note that the function A(0) is nearly constant of value 2£/3 for < 0.1 and tend to 
at = 1. That is to say, in the very early universe, there is a very large cosmological 
constant A ~ O(10 -15 ), but at present day, it will be very small since — > 1. Here we 
adapt the natural units: h = c = G = 1. 

Varying action A with respect to g^, we have field equations 

Rnu - g = — ~~T (4) 

where is the total stress tensor T^ v = i M „ + t^ v . The i M „ and are the stress tensors 
of matter and field respectively, they read 

V = PmQtxv + (Pm + Pm) U^Uy (5) 
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and 



7T0 V Z / 87T 

where p m and p m are energy density and pressure of matter respectively. Varying the 
action A with respect to 0, we have field equation of 0: 

(2u + 3)D0 = 8tt^ + 20 2 A' - 20A - o/0 • - - ^ . (7) 



III. THE COSMOLOGY: BASIC EQUATIONS 



According to the principle of cosmology, = 4>{t). Thus in the Robertson- Walker 
metric 

ds 2 = -dt 2 + a 2 (t) (7^3 + r 2 d9 2 + r 2 sin 2 ed<p 2 ^ , (8) 
there results the following equations: 

77-2 = (a\ 2 _ 8n H^L^A (q\ 

W 3 + 60 2 a0 + 30 Pm 201-0 a 2 ' lJ 

a A c<j0 2 a 87r T 2 

a~3~302~2^0~20~ 30 Pm + 20 1 - ' ^ ^ 

and 

7 ^ , „ ; , 2 , 87r(3p m - Put) r 00 2 2r 00 

+ 3ii0 + — — + 



2(1-0) 2u + 3 2cj + 3(1-0) 2 2^ + 31-0 

= 4/30(1 -0 2 ). (11) 

The energy density and the pressure of the field [] are given by 



1 With the following expressions of and p^, the two equations (9) and (10) can be expressed 
as H 2 = ||(p m + P(j) ) - and | = -|^(p m + P</> + 3p m + 3^) respectively. 
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If cod) 2 ■ •■ 1 Yd) 2 \ 

The energy equation can be found from the identity {R^ u — ^g^R).^ = 0, it is 

p + 3H(p + p) - = 0, {p = Pm + P<t>i P = Pm+P<f>)- (14) 

In turn, the energy equation of matter is given by 

• i. ^ r ( M 1 ft < otj ft \ 

p m + 3i?( Pw + p m ) = - ^_ ^—^Stf— j. (15) 

The terms in the right-hand side are responsible for the creation of matter by gravitation 
characterized by the parameter T. Heating is mainly in the later period of inflation when 
(f) and 4> get larger. When the <f) oscillates after inflation, <p and dj have opposite phases, 
particle creation is insignificant. The value of T is determined by observations of the 
present universe to be: T/£ = 8.482 x 1CT 3 (see section |VI]). The temperature at the end 
of inflation is ~ 10 14 GeV, there is no need of reheating. 

From equation ([0]) we see that the scalar dj moves in a double-well potential 

V{4>) = 0(1 - 2 ) 2 (16) 

which ensures an asymptotic value <j) a = 1 of <p. We note that this double-well is provided 
by the two functions A(0) and uj(dp) together Q . 



IV. THE VERY EARLY UNIVERSE: INFLATION 

A. The creation of the universe 

The quantum creation of the universe starts at t = when the gravitational scalar 
<j)(t) = 0(0) = 0. In fact, the probability for creating the universe is given by || 

P^e- 3< ^ 8y W, (17) 



2 It is given by V'(<f>) = (2c/> 2 A' - 20A)/(2o; + 3). 
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where V is the scalar potential, so, this occurs to the highest local maximum which 
locates at (f> = in the double-well potential. Or, we may argue that || the creation of 
the universe occurs at the zero value of the action (zero action principle), and this is the 
case <fi — d(j) — for the action ([J). 

We are also convinced that the universe is created by gravity from " nothing" . Before 
the creation of universe, there is no any matter but the field 0. If the field takes negative 
value and is trapped in the left well of the potential, the equation d 2 <p/ dr 2 +4/3(1 — (p 2 ) = 
has the instanton solution 

0(r) = tanh y^r, (18) 

where r is the Euclidian time. 



the Euclidean equation 

(l) 2 + ^ 2 -^ < 19 > 



a(r) = - cos x T , (20) 

This solution describes the creation of a closed universe 
1/x- By that time, equation (Q) has the classical solution 

a(t) = — cosh xt ) (21) 

X 



At <f) = 4> = 0, for k = +1, 

has the instanton solution 

here r is the Euclidean time, 
from tunneling. At r = 0, a = 

where 



X = v/2^/3. 

B. The initial condition 



(22) 



The universe after its creation experiences a quantum era whence matter is created. As 
this quantum era is known very little, we shall set the initial conditions at the Planckian 
epoch to such that classical description is adequate hence forward. At the Planckian 
epoch, we can adopt the Planck relations: Ipi = \J~G and t Pi = \fG, and replace the 
gravitational constant G by (f) . It is reasonable to let ao = a(to) — X ■> then the initial 
condition follows as 
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a ~ X ~\ 0o^X- (23) 



The value of p m0 is set to satisfy the causal constraint a H(t ) < 1. In general, p m0 and 
P^q are of the same order. By the way, the exact value of initial p m is of no importance, 
since it is damped away by inflation. 



C. Inflation 



In the very early universe, ~ 0, we see from eqs. C[T2|JT~3D that p^ = — p^. At the 
time t , P<j> an d p m are of the same order, there exists 3p + p < 0, therefore inflation 
commences already. From this time on, p m decreases rapidly by redshift, p^ increases 
and soon dominates, there exists p = —p, then the universe inflates exponentially with a 
Hubble parameter: 

H = X e ■ (24) 

The constant \ e is somewhat different from \ given in (18). Since, in the slow-rolling 
approximation, the term —H<fi/(f) in equation (^) is just —4/3(1 — </> 2 )/3, so that the 2^(3 in 
X((f)) should reduce to (2£ — 4)/3, we set it to be (2£ — 3)f3 to take account of the o;0 2 /60 2 
term. Thus we take 

Xe = v^e/VS , tie = £ - 1-6 . (25) 
The solutions of eqs. (p| JTT|) are 

a(t) = a e Xe(t - t0 \ (26) 



<j>(t) = o e D(t - t()) (27) 

respectively, where 

D = (^1 + 8/(3^-1)^3^/2. (28) 

In the above solutions, we have set to = 0. In equation fl2"?D , we have deleted the decaying 
solution. Substituting the solution ( p7D into equation (pT5|) and integrating, we have 



Tr 2 (3H + rj 2 



4:7i(4:H + 2r) 
7 



<P(t), (29) 



we have deleted the decaying solution also. p m increases with 2 steadily, there is no 
reheating problem. 

When (p(t) rolls down the potential hill, all of <ft, <fi and <ft increase, the energy density 
of the field <fi dominates continuously. When gets larger than 0.4 or so, p^+p^ becomes 
positive and increases as <fi increases. When <fi gets to 0.9 or so, it happens p^ + Sp^ = 0, it 
is the same that p + 3p = 0, then inflation ends naturally. The e-fold number of inflation 
is thus 

N = ^(ln 0.9 - In (f) ). (30) 

We may compare calculations from the above formulae with computer solutions of the 
simultaneous equations (|| [□], |l^) given in Table 1. Here we set 

£ = 7.5 , p = 1.7 x 10" 16 , T = 0.06362. 

The initial conditions according to the above subsection are 

4> = 2.55 x lO" 15 , j> = 2.5 x 10~ 23 , p m = 10" 30 , a = 2 x 10 7 . 
The calculated Hubble parameter, the ratio (/>/</> and the e-fold number are 

Xe = 2.5859 x 10~ 8 , D = 7.9508 x 10~ 9 , N = 109. 



These values calculated from equations (|25|-|30]) can be compared with the numerical 
solutions which are compiled in Table I. 

o The calculated e-fold number is 109 and that from numerical solutions is 110. 

o The Hubble parameter H is nearly constant from Hubble time 7 to 90 because of 
the pressure and the energy density of field (f> are equal in magnitude and opposite in sign 
and > 10 5 x p m . 

o The value D, i.e., <p/4>, are nearly constant in the period mentioned above, that is 
to say, (ft increases exponentially indeed. 

o The initial p m is damped away by 10 Hubble times. In the period 10 — 90 Hubble 
time, the p m {cai) is less than that of numerical solutions a little because of we have ignored 
the second term in the right bracket of equation ( and 1 — is replaced by 1. p m 
increases with (j) 2 in this period. 
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o The temperature at the end of inflation is 3.5 x 10 13 GeV, this is the highest temper- 
ature of the universe. Therefore GUT is never exact and there is no monopole problem 
at all since monopoles are supposed to be produced when the exact GUT is broken. 



Table 1. Numerical solutions of equations (PJTTl |i5|) . 



tHub 



H 



Pr, 



1 

4 
7 

10 
20 
30 
40 
50 
60 
70 
80 
90 
100 

105 
110 



1.9429- lO" 8 
2.4903- 10- 8 
2.5848- 10~ 8 



2.5848 • 
2.5848 • 
2.5848 • 
2.5848 • 

2.5848 ■ 

2.5849 ■ 
2.5850- 
2.5872- 
2.6143- 
2.7722 • 



io- 8 
io- 8 
io- 8 
io- 8 
io- 8 
io- 8 
io- 8 
io- 8 
io- 8 
io- 8 



2.7483 • 10~ 8 
1.1340- IO" 8 



4.9047- IO" 15 
7.4088 • 10~ 15 
3.5084- 10~ 14 

8.8348- 10^ 14 
1.9166- IO" 12 
4.1549- 10~ n 
9.0219- IO" 10 

1.9559 • 10~ 8 
4.2402 ■ IO" 7 

9.207- 10~ 6 
1.9945 ■ 10~ 4 

0.00422587 
0.0753908 

0.76283 
0.939783 



1.1024- IO" 8 
8.4789 • 10~ 9 
7.9534 ■ 10~ 9 



7.9534 
7.9534 
7.9534 
7.9534 
7.9534 
7.9333 
7.9530 
7.9477 
7.8820 
7.4583 



3.9893 • 10~ 31 
5.8882 • 10~ 31 
2.7981 ■ 10~ 30 



1Q -9 
1Q -9 

io- 9 
io- 9 
io- 9 
io- 9 

10 -9 
10 -9 
10 -9 

io- 9 



7.0461 
1.5285 
3.3137 
7.1953 
1.5599 
3.3817 
7.3435 
1.5935 
3.4475 
6.9145 



6.9723 • 10~ 9 
2.3104- 10~ 9 



1.7903- IO" 32 
3.3940 • 10~ 34 



6.9693-10 



1Q -30 
10 -28 

io- 27 

10 -26 

io- 24 

1Q -23 
1Q -22 
1Q -20 
1Q -19 

io- 18 



1.7957 
8.4303 
3.9622 
1.8681 
8.7799 
4.1265 
1.9453 
9.1192 
4.0464 
1.2467 



2.4894 • 10~ 17 
1.4136 ■ 10~ 17 



-43 



IO" 45 

1Q -43 

10~ 40 

io- 37 

1Q -35 
10 -32 
1Q -29 

10~ 27 

1Q -24 

io- 21 



1.8678 • IO" 20 
2.9037- IO" 19 



V. THE UNIVERSE AFTER INFLATION 
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A. Basic equations after inflation 



After inflation, let <fi = 1 — a 2 . We consider the evolution for large time, t ^> (3 1 ( r 
0(1(T 35 ) sec), a < 1. Equations © and (0) become 



« 5?V + (| - 2r ) * 2 + f"™ (31) 

\1 - —(1 - a 2 )\ a + 3H& + Apa = 6/3a 3 - 2/3<x 5 - ya(p m - 3p m ) 

ss 0. (32) 

The energy density and pressure of are 

= ^(4e/?o r2 + ^ 2 + 6^ ( r-6r ( T 2 ), (33) 

P4 , = i- U(2 - + (£ - 2)a 2 + 2Faa - ~ , (34) 



and the equation (|i"5|) becomes 

r 



p m + 3H(p m + p m ) = - {ad + 3H& ). (35) 



7T 



Since the energy density p^ dominates and it is nonrelativistic [] , the universe expands 
by power law about 2/3. From equation fl32"D , we see that a oscillates as 

a(t) = f(t) cos(rot), zu 2 = 4/5/(1 - 2 7 ), (36) 

where 



By equation (|32T), f(t) satisfies 



7= ? - (37) 



' = -WT*T) f - (38) 



3 The exact expression is F 2 = §£/3(cr 2 - (1 - a 2 )Zn(l - a 2 )) + 3 ff M) - + 2f ^ + 

8tt p m _ 2rri- 2 _ 1 

3 T^a 1 a*- 

4 When (j>(t) gets near the point 1, the bottom of the potential well, <\> behaves as particles with 
effective mass ~ lyfft larger than its kinetic energy. 



10 



B. The radiation era: p m = hp 



3 Pm 



By numerical solutions on computer, we see that / ~ j. Let 



and substituting into (j58|),we get 



/(*) = \ , (39) 



Then, substituting f(t) and H(t) into equation ( p5[) and integrating, we get solutions: 



The densities are 



p = (1 ~ 27)7 (42) 

and 

(l_2 7 )(l-8 7 ) 

In the above solutions, the constant c has been determined by equation (|3T|) . 



(43) 



C. The dust era: p m = 

After the time t e<? , the time of radiation-dust equality, p m = 0, the numerical solutions 
show that 

*-s- (44) 



From (pq) , we have 



/(*) ~ t^y) ■ (45) 

After integrating ( P5| ) and determining the integrating constant, we have solutions: 

/1 _ rZ 1 +27/(1-27) 

a(t) = Ji 7^^== T /n „ > cos(rot) ; (46) 
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and 

-47/(1-27) 

(4£ 



1 1 - 87 / t 
6^ 1 - 2 7 \t^ / 



The function a; increases as ~ t 2 . For the given model, for example, at t — lsec, uj 
O(10 73 ). The scalar-tensor gravity agrees with the general relativity exactly. 



VI. THE MODEL 

Now we determine the three parameters (£, /3, T) in the model from observation data. 
The parameter T can be determined by the values of h, T and Q m at present time to as 
follows. 

A. Determining T/£ 

We assume Q = 1 and set h = 0.6, then t = 6.361 x 10 60 = 10.87 Gy and p c = 
1.312 x 10 -123 = 6.765 x 10~ 30 g/cm 3 . The temperature is determined to be T = 1.925 x 
10" 32 = 2.728 fsT and p 7 = 9.039 x 10~ 128 = 4.662 x 10" 34 g/cm 3 . The density of matter 
is assumed to be Q m = 0.4 and so p m = 5.248 x 10~ 124 = 2.706 x 10~ 30 #/cm 3 . Thus the 
temperature of dust-radiation equality is calculated to be 

T eq = 6.651 x 10~ 29 = 9423 K. (49) 



We have the first equation to determine the parameter V and the time t 



eq- 



7(1-27) I — ) = 2.747 x 10 9 , (50) 



where 7 = T/£. On the other hand, relating the expression (|47D for p m (t Q ) to the value 

p m = 1.822 x 10~ 124 gives the second equation: 

-47/(1-27) 



I-87 to_ 

1-27 {teq, 



0.4 . (51) 



Solving the above two equations simultaneously gives 

t /t eg = 573900, t e g = 1.108 x 10 55 = 18940 yr, 

(52) 

7 = 0.008482. 

It should be noted that this result depends only on the adapted values of h,T ,fl m 
but not others. 
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B. The spectra of density perturbations and gravitational waves 



Since the energy density of the field dominates in the evolution of the universe, we 
consider its adiabatic density perturbations. We perform a conformal transformation: 

9iXV = A 2 (4>)~g^, A 2 (0) = i , (53) 

and 



then 



r V 16tt 1 + vT^? # V 16tt ' v ; 



.A = / d 4 x^p g [-^R{g) - \g^V - 2V{$)\, (55) 



where 



The field is a canonical field, then we can calculate the slow-roll parameters 



2 



and 



= C (57) 



77 = — = - . (bl 

8tt V {<t>) £ 



The spectrum of adiabatic density perturbation is 



t2 ,„ 512iry 3 (0) „„„„„ r « 2 



^ W "^ = °' 03395 T- (59) 

In these results, we have noticed that <C 1 for interested k in the range \jhr x Mpc — 
l/SOOOh^Mpc. 

The spectral index is 

dlnd 2 H (k) D . . 

where p = D/if = | (1 + S/^e) 1 / 2 — 1 which depends on £ only. We note that rf/n/c ~ 
Hdt, since at horizon crossing k = aH and a ~ e Ht , H is constant. 
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Inflation also generates gravitational waves, whose relative contribution to the mean- 
squared low multipoles of the CMB anisotropy is 

«(v2*y =10e . (61) 



•57T 



V v{$) ) 



This more correct factor 10 was given by Polarski and Starobinsky 0. For our model, 
r ts = 1.3, this rather large value should be tested by future observations. 

C. COBE 

According to COBE DMR observations of large-angular-scale CMB anisotropy || , we 
normalize the density perturbation spectrum using the result |§: 

Let n s = 0.7f\ i.e., p = 0.3, then 

£ = 7.5 , (63) 

and eqs.(|59|, |61| , |62|) determine 

(3=1.7x 10~ 16 . (64) 
In this model, the energy scale of inflation is 

V l/i = (6.3 x 10 w GeV)e 1/ \ (65) 

To convince the viability of this model, we calculate ag and the large-scale streaming 
velocities cr„(r) which are the three dimensional velocity dispersions of galaxies within 
sphere of radius rh~ x Mpc after the data have been smoothed with a Gaussian filter of 
12 h^Mpc. The results are given Q in Table 2. 



5 Smaller n s gives smaller (£,/?) but larger large-scale streaming velocities. For example, n s = 
0.55, gives (f,/3) =(5.5, 1.7 x 10~ 18 ), a$ = 0.68, ^(40) = 451 km/s, a v (60) = AUkm/s. 

6 The observational data g(j are cr,,(40) = 388(1 ± 0.17) km/s and <r„(60) = 327(1 ± 
0.25) km/s. 
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Table 2. Some parameters of the model 



n s 


^8 


<x,(40) 


<7„(60) 


0.7 


0.69 


357 km/ s 


309 km/s 



We note that in calculating the streaming velocities, we have adopted the CDM tran- 



sition function 11 



T(q) 



In (1 + 2.34g) 
2.34g 



1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71) 4 



(66) 



with q = k/hT s (Mpc~ v ) , and T s = Q h ,exp(—Q B — J -^Q B /Q ) . This manifests itself 
that the gravitational scalar field is just the cold dark matter men seek for, of course, 
there may be other dark matter contained in Q m (in this model, Q m = 0.4, Q B h 2 = 0.019) 



VII. THE PRIMORDIAL NUCLEOSYNTHESIS 



Now, we consider the constraint from primordial nucleosynthesis on scalar-tensor the- 
ories of gravity. Nucleosynthesis is sensitive to the speed-up factor £ n which cannot differ 
from unity by more than ~ 15% without over- or under producing 4 He [12]. The £ n has 
been derived to be 

A(ifn) 



tin 



(67) 



v/l + a% ) 

with ip n being the value of the scalar field at the time of nucleosynthesis and (fo its value 
today. A((p) = I/a/0, cn((p) is defined by [12] 



a(<f) = In A(<p), 



(68) 



a((p) = 



da(ip) 

w 



Here, <p = 
A- 

Therefore 



7T0, since [12] 
1 



d xJ—g 



16n 



R(~g) - 2&^' V - w{<p) 



«M = -^v /r 



(69) 



U{ip) 



A(0) 



(70) 



(71) 
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At the time of nucleosynthesis, A(ip n ) = 1 + O(10 72 ), and at t , a 2 ((p ) is of order 
O(10- 104 ), then 

e„ = l + O(10- 72 ). (72) 
We see that it makes no significant differences with the standard cosmology. 
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